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1979—2017 4 CRU, ERAS、CMED 格 点 降水 
数据 在 青藏 高 原 适 用 性 评估 
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] 1979—2017 年 青藏 高 原 131 个 气象 台 站 实测 降水 资料 ,分 别 从 年 ` 季 尺度 对 CRU ERAS 和 CMFD 3 
种 再 分 析 降 水 资料 在 青藏 高 原 的 适用 性 进行 了 评估 。 结 果 表 明 :(1) CRU ERAS .CMFD 3 种 数据 对 青藏 高 原 年 降 
水 的 模拟 能 力 都 很 强 , 与 观测 值 的 相关 系数 均 超过 了 0.9, 较 观测 降水 量 值 均 偏 大 ;CRU 和 CMFD 的 春季 降水 较为 接 
近 观 测 值 ,CMFD 夏 、 秋 季 降 水 与 观测 值 最 接近 ,3 种 数据 对 冬季 降水 模拟 能 力 均 较 弱 。(2) 从 降水 量 分 布 来 看 ， 


CMED 数据 对 青藏 高 原 年 ,春季 、 夏 季 冬季 降水 的 分 布 模拟 能 力 最 好 ,3 种 数据 对 青藏 高 原 秋季 降水 的 模拟 能 力 具 
有 区 域 差异 性 ,CRU `.CMFD 在 青藏 高 原 西部 的 降水 与 观测 值 较为 接近 。(3) 从 变化 趋势 来 看 ,青藏 高 原 年 春季 IL 
FE 秋季 降水 均 呈 增加 趋势 ,其 中 夏季 增加 幅度 较 大 ,冬季 降水 整体 呈 减 少 趋势 。(4) CRU 数据 对 青藏 高 原 年 春季 、 
夏季 秋季 降水 的 变化 趋势 与 观测 值 较为 一 致 ,其 次 为 ERA5,ERA5 冬 季 降 水 与 观测 值 较为 一 致 。(5) 从 偏差 分 析 
来 看 ,CMFD 数据 与 观测 值 的 年 . 季 降 水 偏差 最 小 ,最 接近 观测 值 .(6) 3 种 数据 与 站 点 平均 的 年 、 季 降水 的 时 间 变 化 
序列 表明 ,CMFD 的 年 春季、 夏季 秋季 降水 与 观测 值 的 变化 最 为 接近 ,其 次 为 CRU;CMFD 冬 季 降 水 较 观 测 值 最 为 


接近 ,但 相关 系数 未 通过 95% 的 显著 性 检验 。 
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降水 作为 影响 青藏 高 原 地 区 自然 生态 系统 最 
活跃 .最 直接 的 气候 因子 ,对 生态 环境 具有 显著 影 
响 ,降水 异常 往 往 极 大 程度 制约 着 生态 环境 的 发 
展 , 在 青藏 高 原 暖 湿 化 大 背景 下 ,青藏 高 原 降水 整 
体 呈 增多 趋势 ,但 区 域 降水 的 空间 差异 显著 ,影响 
因素 不 尽 一 致 “ ,青藏 高 原 不 同 区 域 是 否 具有 一 至 
的 变化 特征 仍 有 待 研究 ,青藏 高 原 地 形 复杂 ,地 域 
广阔 ,部 分 地 区 观测 站 点 稀少 ,往往 代表 性 差 , 单 站 
的 气候 变化 并 不 能 全 面 ,准确 反映 特定 区 域 的 气候 
状况 。 因 此 ,找到 适用 于 青藏 高 原 的 高 分 辨 紊 降水 
资料 是 研究 青藏 高 原 气候 变化 尤其 站 点 稀少 区 域 
气候 变化 的 基本 保障 ,无 论 在 业务 还 是 科学 研究 中 
都 具有 重要 意义 。 

已 有 相当 数量 的 全 球 降水 资料 ,包括 PREC/L 


(NOAA’s Precipitation Reconstruction over Land) 、 
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CMAP(CPC Merged Analysis of Precipitation) .GPCP 
(Global Precipatation Climatology Project) .CRU (Cli- 
matic Research Unit) , ERA (European Center for Me- 
dium-Range Weather Forecasts ) 和 NCEP (National Cen- 
ter for Environmental Prediction ) 全 球 月 降水 资料 等 ， 
这 些 资料 均 已 得 到 广泛 应 用 。 施 能 等 “首先 对 现 有 
不 同 数 据 进 行 了 详细 介绍 ,并 将 其 运用 到 长 江 中 下 
游 的 季风 降水 中 ,发 现 已 经 创建 的 全 球 降 水 量 资料 
有 较 高 的 精确 度 。 随 着 技术 改进 及 融合 数据 多 样 
化 , 自 勇 等 "将 GPCP 数 据 与 中 国 站 点 观测 数据 进行 
了 对 比 , 姜 贵 祥 等 中 利用 1951 一 2010 年 我 国 东 部 
756 个 站 点 实测 降水 资料 ,比较 分 析 了 PREC、CRU、 
APHRO 和 GPCC 4 种 格 点 降水 资料 在 我 国 东部 夏季 
降水 变 率 的 可 靠 性 和 适用 性 ,发现 就 中 国 东 部 夏季 
降水 变 率 的 EOF 模 态 、 主 周期 而 言 ,GPCC 的 适用 性 
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更 高 .更 可 靠 。 在 此 基础 上 , 王 红 丽 等 "将 PREC/L、 
CMAP、GPCP 及 NCEP2 4 种 再 分 析 降 水 资料 与 台 站 
降水 观测 资料 进行 了 多 年 和 逐年 平均 的 年 、 季 月 
降水 量 距 平 的 时 空 相关 分 析 和 方差 分 析 , 比较 了 再 
分 析 资 料 对 该 区 降水 模拟 能 力 的 差异 ,并 对 我 国 东 
部 、 东 北 、 华 南 、 华 北 等 地 的 降水 数据 适用 性 相继 展 
开 了 评估 ,分 析 了 降水 的 变化 趋势 。 

近年 来 ,许多 学 者 对 西北 .青藏 高 原 区 域 降水 
的 适用 性 开展 了 研究 ,如 Wang 等 "通过 比较 青藏 
高 原 63 个 气象 站 点 的 实测 资料 与 MERRA(The Mod- 
ern EraRetrospective-Analysis for Research and Appli- 
cations) .NCEP/NCAR-1 CFSR (NCEP Climate Fore- 
cast System Reanalysis) , ERA- 40, ERA- Interim 及 
GLDAS (Global Land Data Assimilation System) 等 6 
种 再 分 析 资 料 ,发 现 GLDAS 的 降水 数据 与 观测 数据 
最 接近 ;Wang 等 ”依据 站 点 融合 降水 数据 比较 了 
2007 一 2010 年 拉萨 河流 域 4 种 基于 卫星 反 演 的 降水 
数据 (CMORPH 实时 降水 数据 .CMOPRH 站 点 校正 
数据 .TMPA 实时 降水 数据 和 TMPA 站 点 校正 数据 ) 
质量 ,结果 表明 : CMOPRH (CPC MORPHing Tech- 
nique )2 套数 据 的 误差 均 小 于 对 应 的 TMPA 降水 数 
据 , 而 经 过 校正 的 2 套数 据 误差 都 比 二 者 的 实时 降 
水 数据 误差 小 "”。 黄 淆 等 “以 雅鲁藏布江 流域 为 
研究 区 ,验证 了 中 国 地 面 降 水 网 格 数据 .CRU 和 GL- 
DAS 降水 数据 的 精度 ,并 分 析 了 不 同 源 数 据 降 水 量 
年 际 变化 特征 和 概率 分 布 特性 之 间 的 差异 。 孙 畅 
等 “基于 CPCC 和 GPCP 的 降水 数据 及 ERA-interim 
再 分 析 资 料 ,分 析 了 1979 一 2012 年 冬季 青藏 高 原 西 
侧 地 区 降水 的 基本 特征 及 影响 其 年 际 变 率 的 潜在 
子 , 得 出 青藏 高 原 冬 季 降 水 主要 发 生 在 其 西 侧 地 
区 且 为 全 区 变化 一 致 型 ,降水 所 需 的 水 汽 主 要 来 自 
上 游 地 区 ,从 该 区 域 的 西边 界 输入 。 王 坤 讲 等 “ 利 
用 中 国 区 域 地 面 气 象 要 素数 据 集 (CMFD) 和 GLDAS 
数据 集 研究 了 色 林 错 流域 气候 特征 和 时 空 变化 规 
律 ,结果 表明 :气温 的 升 高 和 降水 的 增加 使 色 林 错 
流域 在 近 儿 十 年 来 呈现 显著 的 暧 湿 气候 背景 。 对 
该 流域 的 湖泊 扩张 机 制 .生态 系统 对 气候 变化 的 响 
应 以 及 物候 变化 等 研究 具有 指示 意义 。 

根据 前 人 的 研究 可 知 ,由 于 数据 来 源 和 处 理 过 
程 的 差异 ,不 同 的 资料 对 不 同 区 域 降水 特征 的 描述 
也 不 尽 相同 ,因此 选取 哪 一 种 全 球 再 分 析 降 水 资料 
能 较 好 地 反映 中 国 区 域 降水 特征 ,是 值得 探究 的 问 
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题 。 针 对 青藏 高 原 的 相关 人 研究 ,多数 集中 在 小 区 域 
的 适用 性 评估 ,多 套 再 分 析 资 料 降水 数据 在 整个 青 
藏 高 原 适 用 性 的 评估 比较 缺乏 。 因 此 ,本 文选 取 
CRU ERA5 .CMFD 3 种 再 分 析 资 料 与 青藏 高 原 降水 
观测 资料 进行 比较 ,目的 在 于 找到 一 种 相对 更 适合 
青藏 高 原 降水 研究 的 全 球 再 分 析 资 料 。 一 方面 为 
站 点 稀少 地 区 气候 变化 研究 提供 资料 来 源 ; 另 一 方 
面 可 以 较 好 的 研究 青藏 高 原 气 候 变 化 的 区 域 差 异 ， 
也 便于 与 气候 模式 的 模拟 结果 进行 比较 ,为 青藏 高 
原 气 候 预 测 提 供 一 定 的 理论 依据 。 


1 数据 和 方法 


1.1 站 点 数据 

本 文 所 用 观测 资料 为 1979 一 2017 年 青藏 高 原 
地 区 131 个 基准 气象 台 站 (图 1) 的 逐 月 降水 数据 。 
1.2 CRU 数 据 

东 英 吉利 大 学 气候 研究 中 心 (CRU ) 的 高 分 辩 率 
格 点 化 数据 集 是 目前 全 球 使 用 最 为 广泛 的 近 地 面 
气候 数据 集 之 一 。 该 资料 首先 选择 1961 一 1990 年 
的 气候 平均 值 计算 整个 时 间 段 各 站 点 的 异常 值 , 通 
过 薄板 样 条 插值 法 插值 异常 值 到 格 点 ,再 县 加 气候 
平均 值 最 终 得 到 格 点 资料 。CRU 数据 集 包 含 平 均 
温度 温度 日 较 差 .降水 、 霜 天 和 露天 频率 .水 汽 压 
和 云 量 等 气候 变量 ,水 平分 辨 紊 为 0.5°x0.5*。 本 文 
使 用 的 CRU TS v4.03(CRU Time series version 4.03 ) 
数据 ,时 间 选 取 1979—2017 4F (https://crudata.uea. 
ac.uk/cru/data/hrg/ ) 。 
1.3 ERAS 数 据 

欧洲 中 期 天 气 预 报 中 心 (ECMWF) 作 为 最 早 开 
展 4D -Var 业务 工作 的 机 构 , 其 再 分 析 资 料 研发 也 得 
到 0 了 较 早 的 开展 ,目前 ,ECMWF 再 分 析 资 料 已 经 
发 展 到 了 第 5 代 , 其 在 1979 年 发 布 了 第 1 代 再 分 析 
资料 FGGE, 在 1994 年 发 布 了 第 2 代 再 分 析 资 料 
ERA-15 ,在 2001 年 发 布 了 第 3 代 再 分 析 资 料 ERA- 
40, 在 2006 年 发 布 了 第 4 代 再 分 析 资料 ERA-Inter- 
im,2016 年 发 布 了 第 5 代 再 分 析 数 据 ERA5 的 预览 
版 数据 ,在 2017 年 ECMWF 使 用 4D - Var 20 a 之 际 正 
式 发 布 了 第 5 代 再 分 析 ERA5 数据 ,包括 2010 年 至 
今 的 数据 ,其 水 平分 辩 率 为 0.2$*x0.2$" ,本 文选 取 
1979 一 2017 年 ERA5 月 数据 开展 评估 工作 (https:// 


www.ecmwf.int/en/forecasts/datasets ) 。 
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图 1 青藏 高 原 研 究 区 气象 站 点 分 布 


Fig. 1 Distribution of meteorological station in the Tibetan Plateau 


14 CMFD 气象 要 素数 据 

中 国 科 学 院 青藏 高 原 研究 所 开发 的 中 国 区 域 
地 面 气象 要 素数 据 集 (CMFD) 中 的 格 点 气温 (2 m), 
降水 、 比 湿 .风速 (10 m) 和 太阳 辐射 资料 ,CMFD 数 
据 集 以 国际 上 现 有 的 Princeton 再 分 析 资 料 .GLDAS 
气象 驱动 资料 . GEWEX-SRB 辐射 资料 和 TRMM 降 
水 资料 为 背景 场 ,融合 了 中 国 气 象 局 常规 气象 观测 
数据 制作 而 成 ,水 平分 辨 率 为 0.1°x0.1° ,数据 整体 
准确 度 较 高 ,已 被 广泛 应 用 于 水 文 和 气候 模拟 
研究 。 
1.5 方法 

利用 线性 相关 系数 来 对 比分 析 再 分 析 数 据 对 
观测 数据 的 模拟 程度 ,利用 相对 偏差 平均 绝对 偏 
差 . 均 方 根 误差 等 比 对 再 分 析 数 据 对 观测 数据 的 偏 
离 程度 。 某 一 变量 观测 与 再 分 析 的 相关 系数 计算 
公式 为 : 


n 


D(a- Do; -9) 
R=— (1) 


-aoc 


相对 偏差 为 : 
a;—0 
Bias = 一 x 100 (2) 
Din 


平均 绝对 偏差 为 : 


Mae= Dla, oi/n (3) 
均 方 根 误差 为 ; 
>[e-ac- 可 
RMSE= | 一 (4) 


IP: n 为 有 效 样本 数 ; a P o 分 别 代表 再 分 析 资 料 


和 观测 资料 。 
泰勒 图 能 够 简明 地 总 结 统计 再 分 析 资 料 与 观 


测 值 之 间 的 相 匹 配 程度 。 基 于 余弦 定理 ,泰勒 图 将 
相关 系数 、 均 方 根 误差 和 标准 差 比 率 3 种 统计 数据 
显示 在 一 张 图 中 ,观测 值 被 定 为 沿 横 坐 标 距离 原点 


1 个 单位 的 点 ,距离 原点 的 环形 长 度 为 再 分 析 资 料 
和 观测 值 的 方差 比率 ,代表 再 分 析 与 观测 值 离散 度 
的 相似 性 ;方位 角 的 余弦 代表 相关 系数 ,再 分 析 距 
离 观 测 点 的 长 度 为 再 分 析 与 观测 之 间 的 均 方 根 误 


差 ,距离 观测 点 的 长 度 越 短 ,表明 再 分 析 资 料 越 接 


近 观 测 值 。 


2 结果 与 分 析 


2.1 不 同 源 数据 对 青藏 高 原 降 水 的 模拟 能 力 


为 了 综合 比较 CRU、ERA5 .CMFD 3 种 数据 与 


站 点 观测 降水 资料 的 相似 程度 ,将 它们 的 标准 差 大 
小 和 空间 相关 系数 置 于 泰勒 图 中 (图 2)。 从 年 、 季 


尺度 来 看 ,3 种 数据 对 青藏 高 原 年 总 降水 整体 模拟 
能 力 都 很 强 , 与 观测 值 的 相关 系数 均 超过 了 0.9, 较 
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图 2 1979 一 2017 年 CRU .ERA5 .CMFD 降 水 资料 泰勒 图 
Fig.2 Taylor diagrams of CRU, ERAS, CMFD precipitation 
from 1979 to 2017 in the Tibetan Plateau 


观测 降水 量 值 均 偏 大 。 季 节 尺 度 上 ,3 种 数据 对 青 
藏 高 原 春 、 夏 .秋季 总 降水 整体 模拟 能 力 同样 很 强 ， 
与 观测 值 相关 系数 在 0.73 以 上 , 均 较 观测 值 偏 大 ， 
CRU 和 CMFD 的 春季 降水 较为 接近 ,ERA5 的 春季 
降水 较 观 测 值 偏 大 1.5 倍 以 上 ,CMFD 的 夏季 降水 与 
观测 值 最 接近 ,其 次 为 CRU,ERA5 夏季 降水 较 观 测 
值 同 样 偏 大 1.5 倍 以 上 ,秋季 降水 与 夏季 降水 模拟 
能 力 接近 ,3 种 数据 对 冬季 降水 模拟 能 力 较 弱 。 


(a) 观测 值 


年 降水 量 /mm 
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2.2 不 同 源 数据 的 年 、 季 降水 分 布 

为 了 更 好 的 做 好 格 点 数据 与 观测 数据 的 对 比 
分 析 , 将 站 点 观测 降水 数据 通过 双 线 性 差 值 方法 插 
值 到 0.25°x0.25° 的 格 点 上 ,得 到 观测 降水 的 空间 分 
布 图 ,从 年 降水 的 观测 值 来 看 ,1979 一 2017 年 青藏 
高 原 降 水 自 东 南 向 西北 递减 分 布 ,青藏 高 原 西部 、 
柴 达 木 贫 地 大 部 分 年 降水 量 在 100 mm 以 下 ,祁连山 
区 .西藏 南部 年 降水 量 在 300~450 mm 之 间 ,青藏 高 
原 东 南部 年 降水 量 在 550 mm 以 上 (图 3a) ;对 比 观 
测 值 ,CMFD 数 据 和 ERA5 数 据 能 较 好 的 表征 出 青藏 
高 原 年 降水 量 “ 西 北 少 、 东 南 多 ”的 空间 分 布 特征 ， 
而 CMFD 数据 对 降水 量 级 的 表征 能 力 最 好 ,ERA5 数 
据 对 年 降水 量 的 量 级 偏 大 ,而 CRU 年 降水 量 呈 “ 北 
少 南 多 ”的 空间 分 布 特征 ,从 局 部 来 看 , 柴 达 木 盆 
地 .青藏 高 原 西 北部 及 东南 部 CRU 年 降水 数据 与 观 
测 值 较 为 一 致 (图 3)。 

青藏 高 原 春 季 降 水 观测 值 显示 ,青藏 高 原 中 西 
部 .北部 春季 降水 量 在 50 mm 以 下 ,东南 部 大 部 分 春 
季 降 水 量 为 60~150 mm ,西藏 东北 部 .中 南部 降水 在 
150 mm 以 上 。 对 比 观 测 值 可 知 ,无 论 是 降水 特征 的 
分 布 , 还 是 量 级 的 大 小 ,CMFD 数据 对 青藏 高 原 春 季 
降水 模拟 效果 最 好 ,ERA5 数据 表征 出 青藏 高 原 春 
季 降 水 在 上 某 达 木 岔 地 .青藏 高 原 西 南部 与 观测 值 接 
近 ,青藏 高 原 西北 部 边缘 .东南 部 .东北 部 的 降水 较 
观测 值 偏 大 ,CRU 降水 在 上 某 达 木 贫 地 和 青藏 高 原 东 


(b) CRU 
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图 3 1979 一 2017 年 青藏 高 原 年 降水 量 观测 值 与 CRU ERAS .CMFD 数据 的 分 布 


Fig.3 Annual precipitatin distribution of observation, CRU, ERAS and CMFD from 1979 to 2017 in the Tibetan Plateau 
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南部 与 观测 值 一 致 ,青藏 高 原 中 西部 春季 降水 较 观 
测 值 偏 大 (图 4)。 

从 观测 值 来 看 ,青藏 高 原 夏季 降水 量 大 部 分 在 
50~500 mm 之 间 ,青藏 高 原 西部 .上 某 达 木 岔 地 降水 量 
在 100 mm 以 下 ,青藏 高 原 南 部 、 东 南部 降水 量 在 
250~500 mm 之 间 , 对 比 CRU、ERA5、CMFD 数据 ， 
CMFD 的 夏季 降水 量 分 布 与 观测 值 最 为 接近 ,其 次 


测 值 整体 偏 大 , 约 偏 大 0.75 倍 ,CRU 数据 在 青藏 高 
原 西 部 .上 某 达 木 盆 地 与 观测 值 较为 一 致 ,CRU 数据 
在 青藏 高 原 南 部 的 夏季 降水 较 观测 值 明 显 偏 大 
(AIS). 

青藏 高 原 秋 季 降 水 观测 值 显示 ,秋季 降水 大 部 
分 在 0~80 mm 之 间 ,青藏 高 原 西部 EIA AN eH 
然 是 降水 量 最 少 的 区 域 ,在 10 mm 以 下 ,青藏 高 原 东 


为 ERA5,ERA5 数据 在 青藏 高 原 东 南部 、 南 部 较 观 


(a) 观测 值 


南部 降水 量 在 40~70 mm 之 间 ,CRU ERA5、CMFD 3 


(b) CRU 
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(d) CMFD 
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图 4 1979 一 2017 年 青藏 高 原 春 季 降 水 量 观测 值 与 CRU ERAS .CMFD 降水 量 的 分 布 
Fig. 4 Spring precipitatin distribution of observation, CRU, ERAS and CMFD from 1979 to 2017 in the Tibetan Plateau 
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图 5 1979 一 2017 年 青藏 高 原 夏季 降水 量 观测 值 与 CRU ERAS .CMFD 降水 量 的 分 布 
Fig. 5 Summer precipitatin distribution of observation, CRU, ERAS and CMFD from 1979 to 2017 in the Tibetan Plateau 
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种 数据 对 东南 部 秋季 降水 的 模拟 值 均 偏 大 ,3 种 数 
据 在 柴 达 木 盆地 的 降水 值 均 与 观测 值 较为 接近 ， 
CRU 、CMFD 在 青藏 高 原 西 部 的 降水 值 与 观测 值 较 
为 接近 (图 6)。 

青藏 高 原 冬 季 降 水 量 少 ,观测 值 显 示 ,青藏 高 
原 的 冬季 年 降水 量 约 在 0~20 mm 之 间 ,降水 量 大 值 
区 主要 在 青藏 高 原 西 南 及 东南 边界 ,对 比 可 见 ， 
CMFD 冬季 降 水 分 布 特征 及 降水 量 级 与 观测 值 最 为 


(a) 观测 值 
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接近 ,ERA5 冬 季 降 水 量 在 青藏 高 原 西 南部 、 柴 达 木 
盆地 与 观测 值 较为 一 致 ,其 余地 区 大 部 分 降水 量 较 
观测 值 偏 大 很 多 ,尤其 是 青藏 高 原 西南 部 边界 及 其 
东南 部 ;CRU 冬季 降水 呈 “ 西 多 东 少 ”的 分 布 特征 ， 
这 与 观测 值 明显 不 一 致 ,效果 较 差 (图 7)。 
2.3 不 同 源 降水 数据 的 年 、 季 变化 趋势 

从 年 降水 的 变化 趋势 来 看 (图 8) ,青藏 高 原 大 
部 分 地 区 降水 呈 增 多 趋势 , 约 增多 5~20 mm 
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图 6 1979 一 2017 年 青藏 高 原 秋季 降水 量 观测 值 与 CRU ERAS .CMFD 降水 量 的 分 布 
Fig.6 Autumn precipitatin distribution of observation, CRU, ERAS and CMFD from 1979 to 2017 in the Tibetan Plateau 
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图 7 1979 一 2017 年 青藏 高 原 冬 季 降 水 量 观 测 值 与 CRU ERAS .CMFD 降水 量 的 分 布 
Fig.7 Winter precipitatin distribution of observation, CRU, ERAS and CMFD from 1979 to 2017 in the Tibetan Plateau 
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(a) 观测 值 
K 


年 降水 量变 化 趋势 /[mm.(10a)] 
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二 一 上 
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注 : 打 点 区 域 通过 了 95% 显 著 性 检验 。 下 同 。 
8 1979 一 2017 年 青藏 高 原 年 降水 量 观 测 数 据 与 CRU .ERA5 .CMFD 年 降水 数据 变化 趋势 空间 分 布 
Fig. 8 Annual precipitation variation trend distribution of observation, CRU, ERAS, CMFD from 1979 to 2017 in the Tibetan Plateau 


(10a)…, 其 中 祁连山 区 三江 源 区 ` 改 则 降水 增多 显 
著 ,青藏 高 原 西 南部 .东南 部 边缘 及 柴 达 木 盆 地 西 
北部 呈 减 少 趋势 , 约 减少 5~10 mm“ (10a); CRU K 
水 的 年 变化 趋势 与 观测 值 最 为 接近 ,青藏 高 原 东北 
部 .中 部 .南部 年 降水 呈 显 著 增 多 趋势 , 约 增加 10~ 
30 mm" (10a)! ,青藏 高 原 西 南部 .东南 部 及 北部 的 
柴 达 木 盆地 年 降水 呈 减 少 趋势 ,其 中 青藏 高 原 西 南 
部 年 降水 减少 显著 , 约 减 少 15~25 mm + (10a) ~’; 
ERAS 年 降水 量 整 体 呈 增多 趋势 ,青藏 高 原 南 部 呈 
显著 减少 趋 热 , 约 减少 15 mm: (10a) WE, AK 
部 分 呈 显 著 增 多 趋势 ,尤其 青藏 高 原 西部 大 部 分 地 
区 增多 30 mm (10a)"' 以 上 ;CMFD 年 降水 数据 整体 
呈 增 多 趋势 ,大 部 分 地 区 通过 95% 的 显著 性 检验 ， 


柴 达 木 贫 地 呈 减 少 趋势 外 ,大 部 分 地 区 增多 趋势 通 
过 95% 显 著 性 检验 , 自 北向 南 降 水 增多 逐步 加 强 ， 
约 增多 5~20 mm: (10a)'; ERAS 春季 降水 显示 ,青藏 
高 原 西部 边缘 、 北 部 及 南部 察 隅 地 区 呈 减 少 趋势 ， 
其 中 青藏 高 原 西 部 边缘 降水 减少 趋势 通过 95% 显 
著 性 检验 ,其余 大 部 分 降水 呈 增 加 趋势 ;CMFD 春季 
降水 显示 ,青藏 高 原 春季 降水 呈 显著 增多 趋势 , 增 
多 率 在 3~18 mm*(10a)"! 之 间 , 南 部 局 地 波 密 降 水 减 
少 显 著 , 减 少 速率 在 10~20 mm: (10a). HEIO HE 
较 可 得 ,CRU 春季 降水 整体 呈 增 加 趋势 ,与 观测 值 
较为 一 致 ,其 次 为 ERA5。 

青藏 高 原 东 北部 .上 某 达 木 盆地 .东南 部 青 南 牧 
区 及 西藏 波 密 地 区 夏季 降水 时 减少 趋势 , 约 减 少 2~ 


青藏 高 原 西 南部 琶 拉 木 - 帕 里 、 东 南部 边缘 、 东 部 局 
地 及 南部 波 密 一 带 降水 呈 显 著 减 少 趋势 ,最 大 值 约 
减少 50 mm (10a) "左右, 其余 大 部 分 降水 呈 显 著 增 
多 趋势 ,尤其 是 青藏 高 原 西 北 内 陆 增 多 显著 , 约 增 
多 30~50 mm'(10a)” ,其 次 为 青藏 高 原 东 北部 ,总 体 
而 言 ,CRU 年 降水 变化 趋势 与 观测 值 最 为 接近 。 
青藏 高 原 春季 降水 呈 整 体 增多 趋势 ,其 中 青藏 
高 原 西 南部 普兰 .阿里 地 区 呈 显 著 减 少 趋势 , 约 减 
少 6~12 mm'(10a) ,中 部 及 东部 呈 显 著 增 加 趋势 ， 
约 增加 3~12 mm" (10a); CRU 春季 降水 整体 呈 增 加 
趋势 ,与 观测 值 较为 一 致 , 除 青藏 高 原 西 南部 .北部 


20 mm: (10a) ,其 余 大 部 分 地 区 降水 呈 增 加 趋势 ， 
其 中 青海 省 北部 西藏 改 则 地 区 夏季 降水 增加 显 
著 , 约 增加 10~25 mm: (10a); CRU 夏季 降水 的 变化 
趋势 与 观测 值 较为 一 致 ,青藏 高 原 东 部 红 原 、 阁 尔 
六 一 带 夏 季 降 水 显著 减少 , 约 减 少 10~20 mm + 
(10a), FEAR .由 里 一 带 降 水 显著 增多 , 约 为 23 mms 
(10a)-…;ERA5S 夏季 降水 的 变化 趋势 与 观测 值 的 相 
似 度 仅 次 于 CRU 数据 ,ERA5 夏季 降水 在 青藏 高 原 
西部 增加 20 mm: (10a)"' 以 上 ;从 CMFD 夏季 降水 来 
看 ,青藏 高 原 东南 部 降水 减少 最 为 显著 ,大 部 分 减 
少 速率 在 10~30 mm + (10a)… ,局 部 地 区 减少 速率 在 
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图 9 1979 一 2017 年 青藏 高 原 春季 降水 量 观测 数据 与 CRU ERAS .CMFD 降水 数据 变化 趋势 空间 分 布 
Fig.9 Spring precipitation variation trend distribution of observation, CRU, ERAS, CMFD from 1979 to 2017 in the Tibetan Plateau 


50 mm: (10a) 以 上 ,青藏 高 原 中 西部 和 其 他 数据 一 ”部 呈 减 少 趋势 , 约 减少 3~18 mm'(10a)… ,青藏 高 原 
样 ,为 显著 增加 区 域 (图 10)。 东部 秋季 降水 增多 显著 , 均 通 过 了 95% 显 著 性 检 

青藏 高 原 秋季 降水 趋势 为 南部 及 中 西部 呈 减 ” 验 , 约 增多 3~9 mm'(10a) ,与 观测 值 变化 较为 一 
少 区 域 ( 图 11), 约 减少 3~15 mm: (10a)', 其 余 大 部 。 致 ;ERA5 对 青藏 高 原 秋 季 降 水 变化 趋势 的 模拟 仅 
分 呈 增 加 趋势 ,青藏 高 原 东 南部 的 青 南 牧区 B ”次 于 CRU ,青藏 高 原 北部 降水 增加 显著 , 均 通过 了 
及 红 原 为 秋季 降水 显著 增加 区 域 , 约 增 加 3~9 mm ”95% 的 显著 性 检验 , 约 增加 3~9 mm'.(10a)- ,降水 减 
(10a)'; CRU 秋季 降 水 趋势 为 ,青藏 高 原 南部 .中 西 少 显著 区 域 在 青藏 高 原 南 部 ,尤其 察 隅 一 带 减少 显 


(a) 观测 值 (b) CRU 
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图 10 1979 一 2017 年 青藏 高 原 夏 季 观测 数据 与 CRU ERAS .CMFD 降水 数据 变化 趋势 空间 分 布 
Fig. 10 Summer precipitation variation trend distribution of observation, CRU, ERA5, CMFD from 1979 to 2017 in the Tibetan Plateau 
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图 11 1979 一 2017 年 青藏 高 原 秋季 观测 数据 与 CRU ERAS .CMFD 降水 数据 变化 趋势 空间 分 布 
Fig. 11 Autumn precipitation variation trend distribution of observation, CRU, ERA5, CMFD from 1979 to 2017 in the Tibetan Plateau 


著 , 为 12 mm: (10a)! E ; CMFD 秋季 降水 整体 呈 
增多 趋势 ,青藏 高 原 北 部 降水 增加 明显 , 均 通 过 了 
95% 的 显著 性 检验 ,青藏 高 原 北 部 大 部 分 地 区 降水 
增多 ,在 6~18 mm: (10a)"', 降 水 减少 显著 区 在 青藏 
高 原 东 南端 ,为 15 mm: (10a) U Eo 

观测 值 表 明 ,青藏 高 原 东 北部 、 西 北角 及 东南 
局 部 冬季 降水 呈 增 加 趋势 ,增加 趋势 在 0~2 mm - 
(10a)-… ,青藏 高 原 南 部 .西北 部 主体 冬季 降水 呈 减 
少 趋 势 , 其 中 青藏 高 原 南 部 车 拉 木 、 察 隅 减少 趋势 
明显 ,在 2~10 mm'(10a) ,但 未 通过 95% 的 显著 性 
检验 ;CRU 冬 季 降 水 整体 呈 减 少 趋势 , 除 青藏 高 原 
西北 部 及 东北 部 边缘 冬季 降水 呈 增 多 趋势 外 ,其余 
大 部 分 呈 显 著 减 少 趋势 ,青藏 高 原 西 南部 .南部 边 
缘 冬 季 降 水 约 减少 4~10 mm (10a)-!';ERA5 冬季 降 
水 整体 同样 为 减少 趋势 ,显著 减少 区 域 同样 在 青藏 
高 原 南部 察 阳 一带, 与 观测 值 较为 一 致 ;:CMFD 冬 季 
降水 整体 呈 显 著 增 加 趋势 , 约 2~8 mm (10a), 
减少 区 域 同样 在 高 原 南 部 边缘 (图 12)。 
2.4 不 同 源 年 季 降 水 数据 的 偏差 分 析 

由 于 3 种 再 分 析 降 水 数据 的 水 平分 辩 率 不 一 
致 ,为 了 更 好 的 进行 一 致 性 分 析 , 先 将 观测 值 .CRU、 
ERAS .CMFD 降水 数据 统一 处 理 成 为 0.25°x0.25° 的 
格 点 数据 ,用 再 分 析 数 据 减 去 观测 值得 到 偏差 空间 
分 布 图 ;并 用 青藏 高 原 范 围 内 的 格 点 平均 计算 了 降 


水 整体 的 偏差 , 均 方 根 误差 。 

由 图 13 可 以 看 出 ,年 降水 ,CRU 年 数据 与 观测 
值 相 比 ,降水 量 整 体 偏 少 70.73 mm, 较 观测 值 偏 少 
14.5%, 均 方 根 误差 为 27.66 mm( 表 1) ,从 空间 分 布 
来 看 ,青藏 高 原 西南 部 大 部 分 地 区 偏 多 160~400 mm 
左右 ,其 余 大 部 分 地 区 偏 少 80~240 mm,ERA5 年 降 
水 量 较 观测 值 整体 偏 多 193.11 mm, 偏 多 39.5%, 均 
方 根 误差 为 36.8 mm( 表 1), 除 柴 达 木 伟 地 略 偏 少 
40 mm 左右 外 ,其 余 大 部 分 区 域 均 偏 多 ,尤其 青藏 高 
原 东 南部 偏 多 360 mm 以 上 ,CMFD 年 降水 较 观 测 值 
整体 偏 少 112.55 mm, 偏 少 23.02% , 均 方 根 误差 为 
54.7 mm,3 种 数据 的 年 际 降水 偏差 均 未 通过 95% 的 
显著 性 检验 ;春季 ,CRU 春季 降水 较 观 测 值 整体 偏 
少 10.04 mm(12.33%), 均 方 根 误差 为 12.22 mm, F 
藏 高 原 仍 呈 西南 偏 少 (40~160 mm) .其 余地 区 偏 多 
(0~40 mm ) 的 分 布 态势 ,ERA5 春季 降水 较 观 测 值 整 
体 偏 多 41.51 mm(50.96% ) , 均 方 根 误差 为 7.97 mm 
( 表 1) ,青藏 高 原 大 部 分 地 区 偏 多 0~40 mm, CMFD 
春季 降水 偏差 在 -40~40 mm, 分 布 较为 分 散 , 整 体 较 
观测 值 偏 少 18.66 mm (22.91% ) , 均 方 根 误差 为 
16.32 mm, CRU 和 ERA5 均 通过 95% 的 显著 性 检验 ; 
夏季 ,CRU 夏季 降水 较 观 测 值 整体 偏 少 44.1 mm 
(15.02% ) , 均 方 根 误差 为 21.91 mm, 呈 南 部 偏 多 
(120~400 mm) 北部 偏 少 (80~160 mm) 的 态势 ,其 
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图 12 1979 一 2017 年 青藏 高 原 冬季 观测 数据 与 CRU ERAS .CMFD 降水 数据 变化 趋势 空间 分 布 
Fig. 12 Winter precipitation variation trend distribution of observation, CRU, ERAS, CMFD from 1979 to 2017 in the Tibetan Plateau 


中 柴 达 木 盆 地 略 呈 偏 多 态势 (0~80 mm), ERAS 夏季 
降水 较 观 测 值 整体 偏 多 84.23 mm(28.68% ) , 均 方 根 
误差 为 30.24 mm( 表 1) ,在 青藏 高 原 整 体 偏 多 40~ 
120 mm, CMFD 夏季 降水 较 观 测 值 整体 偏 少 67.98 
mm(23.15%), 均 方 根 误差 为 32.71 mm, 空 间 偏 差 最 
小 ,未 通过 95% 的 显著 性 检验 ,CRU ERA5 均 通过 
95% 的 显著 性 检验 ;秋季 ,CRU 秋季 降水 较 观 测 值 整 
体 偏 少 31.45 mm(30.01%), 均 方 根 误 差 为 8.7 mm, 
青藏 高 原 南部 偏 多 (40~160 mm) .北部 偏 少 (0~40 
mm),ERA5 秋季 降水 在 青藏 高 原 较 观测 值 整体 偏 
多 32.76 mm(31.26% ), 均 方 根 误差 为 8.7 mm, 空 间 
分 布 整体 偏 多 0~40 mm,CMFD 秋季 降水 较 观 测 值 
整体 偏 少 28.92 mm(27.59% ) , 均 方 根 误差 为 12.34 
mm ,空间 分 布 偏差 最 小 ,未 通过 95% 的 显著 性 检验 ， 
CRU 、ERA5 均 通过 95% 的 显著 性 检验 ;冬季 ,CRU 
冬季 降水 较 观 测 值 整体 偏 多 14.63 mm(161.9%), 均 
方 根 误差 为 5.22 mm ,空间 分 布 上 ,整个 青藏 高 原 呈 
西部 偏 多 (0~50 mm) 、 东 部 偏 少 (0~40 mm ) 的 态势 ， 
ERAS 冬季 降水 在 青藏 高 原 整体 偏 多 34.61 mm 
(382.93%), 均 方 根 误差 为 5.6 mm( 表 1), 空 间 分 布 
来 看 ,冬季 降水 整体 较 观 测 值 偏 多 40~120 mm， 
CMFD 冬 季 降 水 较 观 测 值 偏 多 3.13 mm(34.68%), 均 
方 根 误 差 为 5.1 mm, 空 间 偏差 最 小 ,未 通过 95% 的 
显著 性 检验 ,CRU ERAS 降水 通过 了 95% 的 显著 性 


将 格 点 数据 插值 到 131 个 站 点 , 求 取 逐 年 的 站 
点 平均 ,由 图 14 可 见 ,年 降水 、 春 季 降 水 .夏季 降水 、 
秋季 降水 的 观测 值 与 CRU ERA5 及 CMEFD 的 相关 
系数 均 在 0.75~0.98 之 间 , 均 通过 了 95% 的 显著 性 检 
验 , 其 中 CRU、CMFD 与 观测 值 的 变化 很 接近 ,CMFD 
与 观测 值 的 相关 系数 最 好 ,ERA5 的 降水 值 较 观 测 
值 明显 偏 大 ;冬季 降水 ,CRU ERA5 及 CMEFD 的 相关 
系数 分 别 为 -0.03 .0.36、-0.05, 只 有 ERA5 通 过 了 
95% 的 显著 性 检验 ,CMFD 较 观 测 值 最 为 接近 ,但 未 
通过 95% 的 显著 性 检验 ,CRU FI ERAS 较 观 测 值 均 
偏 多 。 


3 讨论 


青藏 高 原 站 点 稀少 ,尤其 是 青藏 高 原 西 北部 几 
乎 没有 观测 台 站 ,青藏 高 原作 为 全 球 气候 变 暖 响应 
的 敏感 区 ,青藏 高 原 的 气候 变化 状况 一 直 都 备 受 关 
注 ,仅仅 依靠 站 点 数据 来 客观 准确、 全面 的 认识 青 
藏 高 原 的 气候 变化 是 不 现实 的 ,寻找 到 在 青藏 高 原 
站 点 密集 .代表 性 好 的 数据 具有 重要 意义 ,也 是 当 
前 气象 业务 发 展 的 必然 要 求 。 关 于 多 源 降水 在 青 
藏 高 原 的 适用 性 评估 工作 ,多 是 围绕 青藏 高 原 局 部 
开展 ,例如 雅鲁藏布江 流域 . 色 林 错 流 域 .青藏 高 原 
西 侧 等 都 是 利用 格 点 数据 揭示 站 点 稀少 地 区 的 气 
候 变化 特征 。 关 于 格 点 降水 数据 在 青藏 高 原 整体 
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图 13 1979—20174F ERAS .CRU .CMFD 年、 季 降水 量 的 偏差 
Fig. 13 Annual and seasonal precipitation difference distribution of ERAS, CRU, CMFD with observation 
from 1979 to 2017 in the Tibetan Plateau 


表 1 1979 一 2017 年 CRU、ERA5、CMFD 年 , 季 降 水 与 观测 值 的 偏差 ,百分比 及 均 方 根 误差 
Tab.1 Difference, percentage and root mean square error of CRU, ERAS and CMFD precipitation compared to 


observation precipitation in annual and seasonals from 1979 to 2017 in the Tibetan Plateau 


与 观测 值 的 降水 偏差 /mm 与 观测 值 的 降水 偏差 百分比 /9% 均 方 根 误差 /mm 
CRU ERAS CMFD CRU ERAS CMFD CRU ERAS CMFD 
春季 —10.04 41.51 -18.66 -12.33 50.96 -22.91 12.22 7.97 16.32 
夏季 -44.10 84.23 -67.98 -15.02 28.68 -23.15 21.91 30.24 32.71 
秋季 -31.45 32.76 -28.92 -30.01 31.26 -27.59 8.70 12.01 12.34 
冬季 14.63 34.61 3.13 161.90 382.93 34.68 5.22 5.60 5.10 
年 -70.73 193.11 -112.55 -14.46 39.49 -23.02 27.66 36.80 54.70 


的 适用 性 评 佑 工作 才刚 刚 开展 ,文中 选取 了 应 用 较 ”能 力 与 谢 欣 汝 等 ”的 研究 结果 一 致 ,CRU 数据 对 
为 广泛 的 3 种 格 点 降水 数据 ,分 析 了 其 在 青藏 高 原 ”青藏 高 原 年 .夏季 冬季 降水 的 模拟 能 力 与 姜 贯 祥 
的 适用 性 ,ERA5 数据 对 青藏 高 原 夏 季 降 水 的 模拟 。 等 的 研究 结果 一 致 ,但 目前 世界 上 的 模式 数据 
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图 14 1979 一 2017 年 观测 数据 与 ERA5 .CRU .CMFD 降水 数据 年 . 季 降 水 量 时 间 变 化 序列 


Fig. 14 Annual and seasonal precipitation timeseries of ERA5, CRU and CMFD from 1979 to 2017 in the Tibetan Plateau 


种 类 较 多 ,评估 结果 具有 一 定 的 相对 性 。 此 外 ,为 
了 便于 与 多 套 格 点 数据 在 空间 上 进行 各 物理 量 的 
对 比 ,对 青藏 高 原 站 点 降水 数据 进行 了 插值 ,因此 
文中 对 青藏 高 原 西 北部 再 分 析 资 料 的 评估 存在 一 
定 的 误差 。 


4 结论 


通过 对 青藏 高 原 年 . 季 站 点 观测 降水 与 3 种 再 
分 析 数 据 的 时 空 特征 、 变 化 趋势 及 偏差 的 比较 分 
析 , 得 出 以 下 结论 : 

(1) 整体 而 言 ,CRU .ERA5 、.CMFD 3 种 数据 对 
青藏 高 原 年 总 降水 量 的 整体 模拟 能 力 都 很 强 ,与 观 
测 值 的 相关 系数 均 超 过 0.9, 较 观测 的 降水 量 均 偏 


大 ;春季 ,CRU 和 CMFD 春季 降水 较为 接近 观测 值 ， 
夏季 ,CMFD 夏季 降水 最 接近 观测 值 ,秋季 降水 与 夏 
季 降 水 模拟 能 力 接近 ,3 种 数据 对 冬季 降水 模拟 能 
力 较 弱 。 

(2) CMFD 数据 对 青藏 高 原 年 春季、 夏季 降水 
特征 分 布 和 量 级 大 小 的 模拟 能 力 都 很 好 ;秋季 ， 
CRU ERAS 、CMFD 对 降水 的 模拟 能 力 具有 区 域 差 
异性 ,3 种 数据 对 青藏 高 原 东 南部 的 降水 的 模拟 均 
偏 大 ,3 种 数据 在 青藏 高 原 北 部 柴 达 木 盆 地 的 降水 
均 与 观测 值 较为 接近 ,CRU .CMFD 在 青藏 高 原 西 部 
的 降水 与 观测 值 较为 接近 ;冬季 ,CMFD 的 降水 分 布 
特征 及 降水 量 级 与 观测 值 最 为 接近 。 

(3) 变化 趋势 :从 年 降水 来 看 ,青藏 高 原 大 部 分 
降水 呈 增 多 趋势 , 约 增多 5~20 mm: (10a) ;春季 降 
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水 整体 为 增多 趋势 ,青藏 高 原 西 南部 的 普兰 .阿里 
地 区 呈 显 著 减 少 趋势 , 约 减少 6~12 mm + (10a)"'; 夏 
Æ ,青藏 高 原 东 部 .东南 部 的 青 南 牧区 及 西藏 的 波 
密约 减少 2~20 mm" (10a)"', 其 余 大 部 分 夏季 降水 呈 
增多 趋势 ;秋季 ,青藏 高 原 南 部 、 中 西部 秋季 降水 约 
减少 3~15 mms (10a)" ,其 余 大 部 分 秋季 降水 呈 增 多 
趋势 ,青藏 高 原 东 南部 的 青 南 牧区 、 玛 曲 及 红 原 显 
著 增 多 3~9 mm: (10a) ;冬季 ,青藏 高 原 东 北部 、 
西北 角 及 东南 部 局 部 冬季 降水 增多 约 0~2 mm， 
(10a) ,青藏 高 原 南 部 .西北 部 冬季 降水 减少 约 2~ 
10 mm’ (10a), CRU 数据 在 青藏 高 原 年 .春季 、 夏 
季 秋季 降水 的 变化 趋势 与 观测 值 较为 一 致 ,其 次 
为 ERA5,ERA5 冬 季 降 水 与 观测 值 较 为 一 致 。 

(4) 偏差 分 析 :CMFD 数 据 与 青藏 高 原 年 、 季 降 
水 观测 值 的 偏差 最 小 ,最 接近 观测 值 , 偏 差 均 未 通 
过 显著 性 检验 。 

(5) 时 间 变 化 序列 :CMFD 数据 与 青藏 高 原 年 
降水 、 春季 降水 .夏季 降水 秋季 降水 的 变化 最 为 接 

近 , 其 次 为 CRU;CMFD 冬 季 降 水 较 观 测 值 最 为 接 

近 , 但 相关 系数 未 通过 95% 的 显著 性 检验 ;ERA5 在 
年 春季 夏季 秋季 及 冬季 的 降水 量 较 观 测 值 均 偏 大 。 
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Assessment of CRU, ERA5, CMFD grid precipitation data for 
the Tibetan Plateau from 1979 to 2017 


WEN Tingting, GUO Yingxiang, DONG Shaorui, DONG Yuanzhen, LAI Xiaoling 
(Qinghai Climate Centre, Xining 810001, Qinghai, China) 


Abstract: In this paper, CRU, ERAS, and CMFD precipitation data for the Qinghai-Xizang Plateau were evaluated 
on annual and seasonal scales and compared with observed precipitation data from 131 weather stations. (1) The 
three data types all showed a strong ability to describe total annual precipitation on the plateau, and their 
correlation coefficients with the observed values were all >0.9. Results from all three models were generally 
higher than observed precipitation. CRU and CMFD results were closest to observed values in spring, while 
CMEFD results were closest to the observed values in summer and fall; all three models showed the weakest 
correlations with observed values for winter precipitation. (2) In terms of precipitation distribution, CMFD data 
had the best ability to simulate plateau precipitation annually, and in spring, summer, and winter. The models 
showed regional differences in describing fall precipitation. CRU and CMFD results were closest to observed 
values on the western plateau. (3) Plateau precipitation increased during spring, summer, and fall, especially in 
the summer, and decreased in the winter. (4) The annual and seasonal variation (spring, summer, and fall) in the 
CRU data was consistent with that of the observed data, and the winter precipitation results from ERAS were 
consistent with observed values. (5) The deviation between CMFD results and observed annual and seasonal 
precipitation was the smallest. (6) Time series analysis showed that CMFD annual, spring, summer, and fall 
precipitation results were closest to the observed values, followed by CRU. CMFD winter precipitation results 
were closest to the observed values, but the correlation coefficient failed the significance test of 95%. 
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